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ABSTRACT 



Context. It is well established that the winds of carbon-rich AGB stars (carbon stars) can be driven by radiation pressure on grains of 
amorphous carbon and collisional transfer of momentum to the gas. This has been demonstrated convincingly by different numerical 
VO , wind models that include time-dependent dust formation. To simplify the treatment of dust opacities, radiative cross sections are 

O^ . usually computed using the assumption that the dust grains are small compared to wavelengths around the stellar flux maximum. 

Considering the typical grain sizes that result from these models, however, the applicability of this small-particle limit (SPL) seems 
, questionable. 

Ph , Aims. We explore grain size effects on wind properties of carbon stars, using a generalized description of radiative cross sections valid 

^0 ' for particles of arbitrary sizes. The purpose of the study is to investigate under which circumstances the SPL may give acceptable 

results, and to quantify the possible errors that may occur when the SPL does not hold. 

Methods. The time-dependent description of grain growth in our detailed radiation-hydrodynamical models gives information about 
dust particle radii in every layer at every instant of time. Theses grain radii are used for computing opacities and determining the 
Q . radiative acceleration of the dust-gas mixture. From the large number of models presented in the first paper of this series (based on 

\—^ • SPL dust opacities; Mattsson et al. 2010) we selected two samples, i.e., a group of models with strong, well-developed outflows that 

TTt. ' are probably representative of the majority of wind-forming models, and another group, close to thresholds in stellar parameter space 

j^ , for dust-driven winds, which are referred to as critical cases. 

Results. We show that in the critical cases the effect of the generalized description of dust opacities can be significant, resulting in 

more intense mass loss and higher wind velocities compared to models using SPL opacities. For well-developed winds, however. 

Cm ' grain size effects on mass-loss rates and wind velocities are found to be small. Both groups of models tend towards lower degrees of 

^ I dust condensation compared to corresponding SPL models, owing to a self-regulating feedback between grain growth and radiative 

acceleration. Consequently, the "dust-loss rates" are lower in the models with the generalized treatment of grain opacities. 

C^^ ' Conclusions. We conclude that our previous results on mass-loss rates obtained with SPL opacities are reliable within a wide re- 

t^^ ' gion of stellar parameter space, except for critical cases close to thresholds of dust-driven outflows where SPL models will tend to 

underestimate the mass loss rates and wind velocities. 

r^ I Key words. Stars: AGB and post-AGB - Stars: atmospheres - Stars: carbon - Stars: circumstellar matter - Stars: evolution - Stars: 

>— ' , mass loss -2 Hydrodynamics - Radiative transfer 
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1 . Introduction lowed by studies of carbon stars including time-dependent (non- 
,,„ , ^ , „ .,,,,,• equilibrium) grain growth (e.g. Fleischer et al 1992; Hofner & 
, Wmds of carbon stars are usually considered to be dust-dnven j^^^^ ^997. ^-^^^^^ ^^ ^j 2OOO) which, despite being based 
^ , wmds. Stellar photons, incident on dust particles, will lead to a ^^ ^^^-^^^^ ^^^^^^^^^ ^U^^^^ ^^ ^^^^^-^^ ^^^-^ 
■ ■ ■ radiative acceleration of the grains away from the star, and, sub- ^-^^ ^^ ^^^^-^ dust-enshrouded carbon stars. In order to ob- 
sequendy, momentum will be transfeiTed to the surrounding gas ^^-^ reasonably realistic results for objects with less optically 
by gas-grain collisions. Pulsation-induced atmospheric shock ^j^^^^ envelopes, however, it is necessary to combine frequency- 
waves contribute significantly to this process by intermittently dependent radiative transfer (including gas and dust opacities) 
creating cool, dense layers of gas well above the photosphere ^^^j^ time-dependent hydrodynamics and non-equilibrium dust 
where dust grains can form and grow efficiently. formation (cf. Hofner et al. 2003). 

Pioneering work on the modelling of AGB star winds was 
done by Wood (1979), focusing on the effects of shock waves. At a point where models of carbon stars are becoming 

and later by Bowen (1988), introducing a parameterized opac- quantitatively comparable to observations as diverse as high- 

ity to describe the dynamical effects of dust formation in the resolution IR spectra (e.g. Nowotny et al. 2010) and spectro- 

circumstellar envelope. These early wind models where fol- interferometric measurements (e.g. Sacutoet al. 2011), it is nec- 

essary to scrutinize a number of underlying physical assump- 

Send offprint requests to: Lais Mattsson tions and approximations. In particular, this concerns the de- 

* e-mail: mattsson@dark-cosmology.dk tailed treatment of dust opacities that are at the core of the 
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Fig. 1. Relevant dust opacity for radiation pressure (combin- 
ing effects of absorption and scattering, see Sect. |2]| as function 
of grain radius and wavelength, computed from refractive in- 
dex data for amorphous carbon by Rouleau & Martin (1991). 
The black contour shows the region where the flux-weighted 
monochromatic opacity exceeds the critical opacity, that is re- 
quired in order for the radiation pressure to balance gravity 
(see Eqs. [12] and [TT] respectively), assuming a Planckian flux 
distribution with Teff - 2700 K and that 30% of the carbon not 
bound in CO condense into carbon dust. 



wind mechanism and also have a direct influence on observable 
properties. A common feature of most detailed dust-driven wind 
models in the literature (including the first paper in this series 
by Mattsson et al. 2010) is that dust opacities are computed un- 
der the assumption that the grain sizes are small compared to the 
relevant wavelengths (defined by the stellar flux distribution), us- 
ing the small-particle limit (SPL) of the Mie theory. In this limit, 
dust opacities are fully determined by the amount of condensed 
material, irrespective of grain sizes, which greatly simplifies the 
modelling because an explicit knowledge of the actual grain size 
distribution in each layer is not required. However, it has been 
shown that grains may grow to sizes where the use of the SPL is 
questionable (e.g., Gail & Sedlmayr 1987, Winters et al. 1994, 
1997, Mattsson et al. 2010). 

The ongoing debate on the mass-loss mechanism of M-type 
AGB stars has recently put the effects of grain size in these ob- 
jects into focus. Using detailed non-grey models, Woitke (2006) 
demonstrated that silicate grains have to be virtually Fe-free in 
the wind acceleration zone, which leads to insufficient radiative 
pressure caused by absorption. Models by Hofner (2008) suggest 
scattering as a possible solution: if conditions in the extended at- 
mosphere allow these grains to grow into the size range of about 
0.1 - lyum, scattering becomes dominant over absorption by sev- 
eral orders of magnitude, opening up the possibility of stellar 
winds driven by scattering on virtually Fe-free silicate grains. 

In carbon stars, on the other hand, the effects of grain size are 
expected to be less dramatic owing to the high absorption cross 
sections of amorphous carbon grains. As shown in Fig.[T] a wide 
range of particle sizes can contribute to driving winds, which 
is different from Fe-free silicates, where small particles are too 
transparent (cf. Fig. 1 in Hofner 2008). For carbon grains with 
radii of about 0.1 - 1/im, however, the SPL of the Mie theory 
may severely underestimate the radiative pressure, with possible 
consequences for the wind properties. The extensive model grid 
presented by Mattsson et al. (2010, hereafter Paper I) shows that 



carbon grains in this size range may be quite common, in partic- 
ular for conditions that allow for efficient grain growth (high C 
abundance, low effective temperatures, slow winds), indicating a 
potential inconsistency with the underlying assumption of SPL 
opacities. 

The main objective of this paper is to establish when the 
small particle approximation can be applied, and to quantify the 
possible errors that may occur in mass-loss rates, wind veloci- 
ties, or dust-to-gas ratios in cases where the dust particles grow 
beyond this regime. For that purpose, we have implemented a 
generalized description of dust opacities in our models, using ac- 
tual mean grain sizes and corresponding radiative cross sections 
that are valid beyond the SPL (see Sect.|2]l. From the model grid 
in Paper I (based on the SPL) we select a subgroup of models 
that show large grains and/or slow winds, and which we there- 
fore can expect to be noticeably affected by the assumptions 
about dust opacities. Re-computing these models with the newly 
implemented, generalized treatment of grain opacities gives esti- 
mates of the errors introduced by using the small particle limit. It 
should be noted, however, that the extreme cases discussed here 
are not necessarily representative of the majority of the wind 
models in Paper I, but that they rather highlight grain size as a 
potentially critical property, and presumably give an upper limit 
of the errors. 

2. Modelling method 

The results discussed in this paper are based on dynamic atmo- 
sphere and wind models that combine non-equilibrium dust for- 
mation and frequency-dependent radiative transfer (taking both 
molecular and dust opacities into account). The effects of stellar 
pulsation are simulated by a variable inner boundary ('piston' 
with accompanying luminosity variation) below the stellar pho- 
tosphere. The general modelling method has been described in 
detail by Hofner et al. (2003) and Mattsson et al. (2007, 2010). 

The new models presented here include a description of dust 
opacities that is applicable to grains of arbitrary size, in contrast 
to our earlier carbon star models, which used a simple limit- case 
of the Mie theory - valid for particles much smaller than the 
relevant wavelengths only - irrespective of the actual emerging 
particle sizes. The more general method of computing radiative 
cross sections can affect both the radiative energy transfer (tem- 
perature structure) and the radiative pressure on dust grains, and, 
consequently, the acceleration of the wind, if dust particles grow 
beyond sizes where the SPL holds. Below, we discuss the newly 
implemented description of dust opacities in detail. 

2.1. Dust opacities: dependence on grain size 

A crucial quantity in the following discussion is the radiation 
pressure efficiency factor Qi-p, defined as the ratio of the cor- 
responding radiative cross-section, Crp, to the geometric cross- 
section of a grain, i.e.. 



2,p(flgi-, A) = 



Crp(agi,A) 



(1) 



assuming spherical grains with radii Ogi. The cross section de- 
termining radiative pressure is a combination of absorption and 
scattering cross-sections (Cabs and Csca, respectively). 



^rp ~ *-^abs "I" \J^ ,?sca/ *-^sca» 



(2) 



with ^sca denoting the mean cosine of the scattering angle, where 
^sca - 1 corresponds to pure forward scattering (see, e.g.. 
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iKriigel 2003h . These quantities can be derived from refractive 
index data of relevant grain materials using the Mie theory. With 
the definitions given above, the opacity that determines the ra- 
diative pressure on an ensemble of dust grains embedded in a 
gas - with p denoting the mass density of the gas-grain mixture 
- can be expressed as 



flgrGrp(flgr, '^) "(flgr) dOg,. 



(3) 



where n(agr) dug^ is the number density of grains in the size in- 
terval t/flgr around Ogi. By defining Q' = 2ip(agr, /l)/flgi and its 
grain-size average 



(g;) 



2m(«gr,^)flgi-n(flgr)^flgr 

Jo 





(4) 



the opacity can be reformulated (without loss of generality) as 



.dust 



P 



{Q'^) 



flgj n(flgi-) da„ 



(5) 



which is a more suitable form for the following discussion. 

In our models the dust particles at distance r from the stellar 
center, at time t, are described in terms of moments Ki(r, t) of the 
grain size distribution function «(flgi-, r, t). 



Kiir, t) oc 



Jo 



"gi- n(flg 



, r, t) da„ 



(/ = 0,1,2,3). 



(6) 



It follows from this definition that A'o is proportional to the to- 
tal number density of grains (the integral of the size distribu- 
tion function over all grain sizes), while Ki, K2, and ^3 are re- 
lated to the average radius, geometric cross-section and volume 
of the grains, respectively. The equations determining the evolu- 
tion of the moments Ki(r, f) (including nucleation, grain growth, 
and evaporation; cf. Gail & Sedlmayr 1988, Gauger et al. 1990) 
are described in detail in previous papers (see Hofner et al. 2003 
and references therein). 

Regarding the computation of dust opacities, the integrals 
over grain size in Eq. (|5]l and in the denominator of Eq. (HJi are 
given by the moment K^, while, in general, the time-dependent 
local grain size distribution in each layer of the model has to be 
known to evaluate the remaining integral in Eq. (HJl, involving 

e;p(«gr,^). 

2.2. The small-particle approximation and its limitations 

In the limit case of particles that are much smaller than the rele- 
vant photon wavelengths, i.e., 27rflgr <sc A, however, the problem 
of computing dust opacities becomes much simpler. According 
to the Mie theory (see, e.g., Bohren & Huffman, 1983), the ab- 
sorption and scattering efliciencies for small grains behave like 
Gabs "^ flgr and 2sca '^ «!• In this limit, absorption dominates 
over scattering, implying that Qrp ~ Gabs, and, consequently, 
that Q' X! 2abs/flgr becomes independent of the grain size, 
making the integration in Eq.© trivial. Therefore, provided that 
27rflgj « A holds for all relevant wavelengths, the opacity can be 
reformulated as 



C^-GlbsW 



I «gi« 

Jo 



(ag,)da„ oc Q' K3, 



(7) 



This expression only depends on the total amount of mate- 
rial condensed into dust (given by ^^3). Consequently, explicit 
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Fig. 2. Histogram of resulting mean grain sizes (derived from 
moment Ki at the outer model boundary; see Eq. ( fTOl ) and text) 
for wind-forming models taken from Paper I, with M* = IMq 
and Aup = 4 km s"', spanning a range of stellar luminosi- 
ties, effective temperatures and carbon abundances (see Table 2 
in Paper I). The vertical dashed line marks a grain radius of 
4 • 10"^ cm where deviations in the opacity from the small- 
particle limit (SPL) may exceed 10% at wavelength A - lyum 
(seeFig.[3]l. 



knowledge of the grain size distribution is not required, which 
greatly simplifies the modelling. Therefore, many models in the 
literature, including those presented in Paper I, have used the 
SPL to describe dust opacities. 

In view of the resulting mean grain sizes for the model grid 
presented in Paper I, however, it is necessary to investigate the 
possible effects of size-dependent grain opacities beyond the 
small particle limit on mass-loss properties of C-type AGB stars. 
A comparison of Fig|2] (showing a histogram of typical grain 
sizes in Paper I) with Figl3] (showing the deviations from small 
particle cross sections as a function of grain size) demonstrates 
that typical grain sizes are in a range where the small-particle ap- 
proximation may lead to a considerable underestimation of the 
radiative pressure. 

2.3. Dust opacities beyond the small-particle limit 

In principle, the size distribution function n(flgr, r, f) can be re- 
constructed from the moments Ki(r, t), allowing for a general 
treatment of dust opacities (Eqs. (|5]l and (HI), but this involves 
a considerable computational effort, well beyond the scope of 
this paper. Instead, we test the influence of size-dependent dust 
opacities on AGB winds with several descriptions that 

- do not require explicit knowledge of the size distribution, 

- are based on the actual moments Ki(r, t), and 

- are not restricted to the small particle limit. 

All these descriptions are in the general form of Eq. (O, which 
means that they are dependent on the total amount of material 
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Fig. 3. Radiative pressure efficiency factor 2rp and its compo- 
nents gab,s and Qscn, relative to the corresponding small-particle 
limit (SPL) values 2spl, as functions of grain radius at A - I fim 
(i.e., near the stellar flux maximum). Data are given for amor- 
phous carbon dust, taken from Rouleau & Martin (1991) and 
the Q's are calculated using the Mie theory for spherical parti- 
cles (programme BHMIE from Bohren & Huffman 1983, modi- 
fied by Draine, www.astro.princeton.edu/draine/scattering.html). 
The vertical dashed line marks a grain radius of 4- 10"^ cm where 
deviations in the opacity from the SPL may exceed 10% at wave- 
length A - I /im. 



actually condensed into grains in each layer at a given instant 
(integral over grain size, proportional to Kj,), 



^f''^ {Q'(r,t))Ki(r,t), 



(8) 



but they use different approximations for {QL(r, t)}. A common 
feature of the new models discussed here is that the size-average 
of Q' (as defined in Eq. |4]i is approximated by the value of Q' 
for an average grain size, agi-{r, t), in each layer, i.e., 

iQ'rp(r,t))^Q'^(a,,(r,t),A). (9) 

More precisely, we consider the following cases: 

(A) Models where a fixed grain radius Ogi is used when com- 
puting Q' (agi-,A) while dust formation (determining Ky,) is 
modelled as usual. We study three cases, where Ogi is taken 
to be either deep in the SPL (10^^ cm), or 'optimised' such 
as to obtain maximum Qrp around the stellar flux maximum 
(3.55 ■ 10"^ cm), or very large (10"-^ cm). In the 'optimised' 
case, 2ip is about four times higher compared to the SPL at 
lyum (see Fig.[3]l. Models with a fixed grain radius are con- 
sidered mainly for reference. 

(B) Models where 5gr is a mean grain radius derived from the 
moments Ki(r, t) (with / = 1,2, 3), i.e.. 



agr(r, t) = (a' } 



i .i/i 



i/i 



where Omon is the monomer radius of the grain material. 
Which of the three cases gives the best approximation of the 
true (Q'} (defined in Eq. |4| is impossible to say without 
prior knowledge of the properties of n(flgr). Hence, we tested 
all three cases. The different mean radii Sgr will be referred 
to as Ki , K2 or Ki, mean, respectively. From a physical point 
of view, the .^1 mean is simply the mean grain radius result- 
ing from the size distribution as such, while the K2 and K3 
means represent grain radii corresponding to the mean grain 
surface and the mean grain volume, respectively. 

In all cases the quantity Qi-pia^rir, t). A) was calculated using 
the Mie theory for spherical particles of arbitrary size (using 
the programme BHMIE from Bohren & Huffman 1983, modi- 
fied by Draine, www.astro.princeton.edu/draine/scattering.html) 
and refractive index data for amorphous carbon dust taken from 
Rouleau & Martin (1991)0 

For completeness sake, we also mention here that the size de- 
pendence of the opacity relevant for determining the grain tem- 
perature, i.e. QihAa„^{r,t),A), is treated in a similar way as the 
radiative pressure. 



3. Definitions and selection of models 

An average grain radius less than 4 • 10"^ cm is referred to 
throughout as "small" regarding radiative cross-sections. This 
value corresponds to the lower limit of the size range where the 
actual 2rp (at wavelength A - \ ^m) may deviate by more than 
10% from the value given by the small particle limit (cf. Fig.O. 
Below we will also refer to dust opacities as being of type A 
or type B:, i.e., those with a fixed grain size (type A; for test- 
ing) and those using varying characteristic grain sizes based on 
mean values computed at each time step and spatial grid point 
throughout the simulation (type B; see Sect.|2]for details). 

We present re-computations of two groups of models se- 
lected from Paper I, using the modified version of our code (see 
Sect. |2]i. The first group of models (numbers 1-12 in Tab. [U, 
referred to as "critical cases" in the following, we expect to be 
significantly affected by including grain-size effects. By "criti- 
cal" we mean that the flux-mean dust opacity is comparable to 
the critical opacity. 



4;rcGM* 



(11) 



which corresponds to a ratio of unity for (outwards-directed) ra- 
diative and (inwards-directed) gravitational acceleration for stel- 
lar luminosity L* and stellar mass M* (c and G denote the speed 
of light and the constant of gravity, respectively). In other words, 
critical cases are defined by a situation where radiation pressure 
on dust is close to the value required for balancing gravity. In 
practice this corresponds to models with slow winds in which 
the dust grains have time to grow larger than usual, and mod- 
els near some mass-loss threshold in stellar parameter space for 
which a slight increase/decrease of Qrp could enable or prevent 
wind formation. 

In addition, a "control group" consisting of 12 models with 
strong, well-developed outflows that also show relatively large 



' Note that for all the descriptions of {Q' ) listed above the SPL of 
the Mie theory is recovered if the assumed or actual grain sizes (in type 
A or type B opacities, respectively) are much smaller than the relevant 
wavelengths (defined by the stellar flux distribution), i.e. Ina^^ « A. 



(10) 



In particular, opacities of type A with 5gj 
comparable with Paper I. 



10 cm should be directly 
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average grain sizes (models 13-24 in Table [T]) was selected and 
recomputed for comparison. These models are producing big 
dust grains according to the definition above (see Fig.|2]i, but the 
effects of grain-size dependent opacities may not be that signif- 
icant, because the momentum transfer efficiency (from radiation 
to dust and gas) in these cases is near the theoretical maximum, 
i.e., the single-scattering limit, which corresponds to the mass- 
loss rate, M ~ Li, Uo^Jc, where L* is the luminosity, Mqui is the 
flow speed over the outer boundary and c is the speed of light. 
The optical depth of the wind is high and the wind speed cannot 
be much affected by an increase of the radiative pressure effi- 
ciency factor Qi-p. 



for models with maximized type A dust opacities compared to 
SPL models (see Fig. |4]i. This is a consequence of dust grains 
having less time to grow; the flow is generally faster and there- 
fore the grains pass through the dust formation zone in a shorter 
time. In general the mean degrees of dust condensation (fc) are 
lower, which makes the "dust-loss rates" several times lower. 

The "control group" of models (13-24) with strong, well- 
developed dust-driven outflows (when using the SPL) are, as ex- 
pected, not much affected by maximizing the dust opacity, as 
far as mass-loss rate and wind velocity are concerned. Those are 
also models that presumably are quite representative of the ma- 
jority of wind-producing models in Paper I. 



4. Results and discussion 

4. 1 . Basic tests and constraints 

To test the modified code, we have tried to replicate the results 
from Paper I by adopting a small fixed grain radius Ogr = 10"^ 
cm, when calculating Q' , which should be well within the small 
particle regime. The average mass-loss rates, wind speeds and 
mean degrees of dust condensation that we obtained are indeed 
almost exactly the same as in Paper I, which indicates that the 
modified code is working properly. 

In the opposite limit, i.e., when the particles are much larger 
than the wavelengths under consideration, Qrp approaches a con- 
stant value (see, e.g.. Fig. 11 in Paper I), and, consequently, 
Q' oc 1 /flgj. From Eq. ^ it can therefore be deduced that for 
a fixed total amount of dust material per volume (represented 
by the integral over grain size), a:JJ"^' k I/agr for 2;ragr » A. 
Because the total amount of grain material per volume is lim- 
ited by the availability (abundances) of the constituting chemi- 
cal elements, there is a limiting maximum grain-size where the 
flux-mean opacity necessarily drops below the critical opacity 
/Ccrit and radiative pressure alone cannot overcome gravity. 

Fig. [1] illustrates the dependence of the dust opacity on both 
grain size and wavelength, and shows an estimate of which grain 
sizes will be relevant for driving winds. The colour scale rep- 
resents the quantity A-jJ"^V'<'crit and the black contour marks the 
region where the monochromatic flux-weighted opacity 



K,-p(A, Qgr) 



7tABAT,b) 



(12) 



exceeds the critical opacity /Cdit, assuming a Planckian flux dis- 
tribution B,i with Teff - 2700 K and that 30% of the carbon 
not bound into CO is condensed into grains (with a free carbon 
abundance ec - eo = 3.3 ■ 10""*, M, = 1 Mq, L* = 7000 Lq and 
cr denoting the Stefan-Boltzmann constant). An upper limit in 
grain sizes relevant for driving a wind is clearly apparent from 
this plot. Adopting a large fixed grain radius (ugr = 10""^ cm) 
should therefore prevent the formation of dust-driven outflows, 
which is confirmed by our detailed models. 

4.2. Recomputed models witli optimized type A opacities 

By "optimizing" the grain radius used in the opacities for max- 
imum Qi-p, we obtain much faster winds and higher mass-loss 
rates for the critical cases (Models 1-12) compared to using the 
SPL (see Tables [T]&|2] and Fig. |4l left column). Four out of five 
cases where SPL models have no resultant wind do indeed have 
a considerable outflow when using an optimized Q^p. which is 
roughly a factor of five higher compared to the SPL value. 

The actual average grain radius derived from the moment 
equations for dust formation, i.e., (flgr) - Omon ^i /Kq, is smaller 



4.3. Recomputed models with type B opacities 

In reality, not all grains can obtain the same radii, and especially 
not the radius that would maximize Q^. The models with type 
B dust opacities take into account that Qi-p changes in time and 
space owing to variations of the grain-size distribution, by using 
a value of {Q'(r, t)) based on different moments of the grain- 
size distribution (cf. Sect. |2]i. For the critical-case models, this 
results in much faster wind speeds and more intense mass loss 
compared to using the SPL (see Fig.|4]i. The "control group" of 
models with strong winds are again much less affected in these 
respects. 

The results are generally quite similar to those of the models 
with maximized type A opacities, which suggests that the aver- 
age grain size tends to be such that 2rp is actually close to max- 
imized in the relevant part of the spectrum. The average grain 
radii are in many cases larger than according to type A mod- 
els, but still smaller than the grain radii inferred from the K\ 
moment of the corresponding SPL models in Paper I. The fact 
that the grain radii tend to be smaller when using the general- 
ized description of dust opacities does not mean that we are ap- 
proaching the small-particle region again (where the SPL holds 
exactly), but reflects the effects of self-regulation in the wind 
mechanism. When dust grains grow beyond the SPL regime, the 
radiative acceleration becomes more efficient and they are then 
likely to move away from the dust formation zone faster, which 
means that they cannot continue to grow. If the momentum trans- 
fer from the radiation field to the dust-gas mixture is not suffi- 
cient to sustain an outflow, the dust grains may on the other hand 
continue to grow, which means that small particles (experiencing 
too little radiation pressure) may grow until they reach optimal 
size. 

The mean degree of dust condensation {fc) is typically much 
lower for critical-case models with type B opacities than for the 
corresponding SPL models of Paper I, but slightly larger than 
for models with type A opacities (see Fig.|4]i. There is one crit- 
ical Kt, mean model (Model 5) that stands out from the rest and 
shows a {fc) value that is several times higher than for the cor- 
responding models with maximized type A opacities, or when 
{QL) is computed with the K\ and Ko mean grain radii (again, 
see Fig. nil. The reason is that this particular model shows no net 
outflow, which means that grain growth is not stopped by falling 
densities as would be the case in a wind. 



4.4. Thresholds for dust-driven winds 

In Paper I we argued that for a realistic description of mass loss it 
is also crucial to know in which parts of stellar parameter space 
dust-driven mass loss cannot be sustained. Thresholds for dust- 
driven outflows originate from the simple fact that there exists 
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Fig. 4. From top to bottom: Mass-loss rates, wind speeds, mean degrees of dust condensation, dust-to-gas ratios and mean grain 
radii for models with type A opacities ("optimized Qip", first column) and type B opacities (using actual grain sizes based on the 
moments K\, K2 and Ky, columns 2, 3 and 4, respectively) vs. the corresponding quantities in SPL models. The dashed lines show 
the case of equal values. 



a critical radiative-to-gravitational acceleration ratio and thus iMattsson et al. 2007bl) . In Fig. |5] a threshold in luminosity is 
critical values of stellar parameters and element abundances for shown as an example. The transition regions in parameter space 
which the flux-mean opacity equals A:crit ( iDominik et al. 1990; from windless models to strong outflows are not adequately cov- 
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Fig. 5. Mass-loss rates (left panel), wind speeds (middle panel) and mean degree of dust condensation near a "threshold" as function 
of luminosity. All other stellar parameters are kept the same (M* = IMq, T^ff = 2600 K, log(C - O) = 8.80). 



ered with the grid spacing chosen in Paper I. To better resolve the 
transition region and to quantify the effect of grain-size depen- 
dent opacities on the threshold, we have computed a number of 
additional models close to the expected mass-loss threshold tak- 
ing smaller steps in log(L*) while keeping all other parameters 
constant. It is evident from the selected models that if the dust 
grains have the most favorable size (in terms of achieving the 
highest possible opacity, i.e., using optimized type A opacities), 
it will be easier to sustain an outflow near a threshold. It is also 
clear that the wind speed near a mass-loss threshold is signifi- 
cantly affected if the typical grain size is in the optimal range. 
In Fig. |5] the wind speed is increased by approximately a factor 
of two relative to the corresponding SPL models. The mass-loss 
rate differs most at the lowest luminosity, log{Li, / Lq) = 3.55, 
in which case grain-size effects mean the difference between no 
outflow and a sustained dust-driven wind. At the upper end of 
the tested luminosities, i.e. log(L*/LQ) - 3.85, the mass-loss 
rate is about 50% higher than in the SPL case, while it is almost 
unaffected for the intermediate values of L*. The mean degree 
of dust condensation /^ is roughly a factor of two lower com- 
pared to the SPL models, like in the critical case models with 
optimized type A opacities discussed above. 

5. Conclusions 

In the first paper of this series we presented a large grid 
of frequency-dependent dynamic models for atmospheres and 
winds of C-type AGB stars, with the main purpose of provid- 
ing a more realistic description of dust-driven mass loss as in- 
put for stellar evolution models. One of the underlying assump- 
tions of the models, i.e., that dust opacities can be described 
with the small-particle limit (SPL) of the Mie theory, however, 
turned out to be questionable with typical emerging grain sizes 
in a range where the SPL may severely underestimate the actual 
grain opacities. 

Introducing a generalized description of radiative cross sec- 
tions valid for arbitrary particle radii, we have explored grain 
size effects on the wind properties of carbon stars. The time- 
dependent description of grain growth used in our models readily 
gives particle radii corresponding to various means of the grain 
size distribution in every layer at every instant of time. To keep 
the computational effort at a level that will allow the construc- 
tion of large model grids, we have used descriptions of the dust 
opacities based on these local mean grain radii. 



From the large number of models presented in Paper I we 
selected two samples, i.e., a group of models with strong, well- 
developed outflows and another group close to thresholds for 
dust-driven winds in stellar parameters space, referred to as 
"critical cases". For the first group, which is presumably rep- 
resentative of most wind-forming models in Paper I, the effects 
of grain size on mass-loss rates and wind velocities are small, 
whereas the critical models show more intense mass loss and 
(in some cases significantly) higher wind velocities when us- 
ing the new, generalized description of dust opacities. Models 
in both groups tend towards lower dust-to-gas ratios, illustrat- 
ing a self-regulating feedback between grain growth and the 
increased opacity (and, consequently, higher radiative acceler- 
ation) per dust mass. Therefore, in general, the "dust-loss rates" 
can be expected to be lower than in Paper I. 

Extrapolating from the samples investigated here, it seems 
that the mass loss rates given in Paper I are reliable, within the 
limits of current theoretical and observational uncertainties, ex- 
cept for models close to thresholds for dust-driven outflows in 
stellar parameter space where mass loss is probably underesti- 
mated. 

Given the results presented here, a full implementation of 
grain-size-dependent opacities seems to be important as future 
work. Furthermore, the actual sizes of dust grains may be of 
great importance for theoretical spectra of dynamic atmosphere 
models of carbon stars, which we plan to study in the near future. 
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Table 2. Same as in Table[T] but using the "2rp-optimized" grain radius in the dust opacities. 
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Table 3. Same as in Table[T] but using the first (Ki) mean grain radius in the dust opacities. 
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2400 


8.80 


524 


I.66E-05 


3.03E+01 


1.23E-0I 


7.00E-04 


I.56E-05 


16 


4.00 


2400 


9.10 


524 


2.02E-05 


4.61E+01 


6.77E-01 


7.69E-03 


7.88E-06 


17 


3.85 


2600 


8.80 


393 


4.15E-06 


2.61E+01 


2.64E-0I 


1.50E-03 


1.24E-05 


18 


3.85 


2600 


9.10 


393 


6.94E-06 


4.64E+01 


5.58E-01 


6.34E-03 


7.06E-06 


19 


4.00 


2600 


8.80 


524 


9.22E-06 


2.97E+01 


3.04E-0I 


I.73E-03 


I.33E-05 


20 


4.00 


2600 


9.10 


524 


I.22E-05 


5.I4E+01 


6.25E-01 


7.10E-03 


7.24E-06 


21 


4.00 


2800 


8.80 


524 


6.32E-06 


2.76E+01 


2.46E-0I 


I.40E-03 


1.22E-05 


22 


4.00 


2800 


9.10 


524 


8.40E-06 


5.28E+01 


5.22E-01 


5.93E-03 


6.30E-06 


23 


4.00 


3000 


8.80 


524 


3.87E-06 


2.08E+01 


2.00E-0I 


1.I4E-03 


I.09E-05 


24 


4.00 


3000 


9.10 


524 


7.76E-06 


4.51E+01 


4.65E-01 


5.28E-03 


6.33E-06 



10 L. Mattsson and S. Hofner: Mass Loss from Carbon Stars as Function of Stellar Parameters. 11. 

Table 4. Same as in Table[T] but using the second {K2) mean grain radius in the dust opacities. 



Mod. 


log(L,) 


Teff 


log(C-O) 


P 


<M> 


<"out> 


<fc> 


iPdlPg) 


(agr) 




[Le] 


[K] 




[days] 


[Moyr-'] 


[kms-'] 






[cm] 


1 


3.70 


2400 


8.20 


295 






4.06E-01 


5.52E-04 


2.32E-04 


2 


3.70 


2400 


8.50 


295 


1.98E-06 


1.71E+01 


l.llE-01 


3.01E-04 


2.00E-05 


3 


3.85 


2400 


8.20 


393 


2.24E-06 


4.97E+00 


1.19E-01 


1.62E-04 


4.31E-05 


4 


3.85 


2400 


8.50 


393 


4.30E-06 


1.72E+01 


1.08E-01 


2.93E-04 


2.00E-05 


5 


4.00 


2400 


8.20 


524 


7.20E-06 


8.12E+00 


6.83E-02 


9.28E-05 


2.61E-05 


6 


4.00 


2400 


8.50 


524 


8.74E-06 


1.59E+01 


9.65E-02 


2.62E-04 


1.73E-05 


7 


3.70 


2600 


8.50 


295 


5.80E-07 


1.45E+01 


7.60E-02 


2.06E-04 


1.83E-05 


8 


3.85 


2600 


8.20 


393 


9.23E-07 


4.97E+00 


l.lOE-01 


1.49E-04 


3.95E-05 


9 


3.85 


2600 


8.50 


393 


1.96E-06 


1.87E+01 


8.41E-02 


2.28E-04 


1.97E-05 


10 


4.00 


2600 


8.20 


524 


3.59E-06 


9.11E+00 


8.13E-02 


l.lOE-04 


3.60E-05 


11 


4.00 


2800 


8.20 


524 


- 


- 


- 


- 


- 


12 


3.70 


3000 


8.80 


295 


- 


- 


- 


- 


- 


13 


3.85 


2400 


8.80 


393 


5.89E-06 


3.02E+01 


2.99E-01 


1.70E-03 


1.33E-05 


14 


3.85 


2400 


9.10 


393 


l.lOE-05 


4.40E+01 


6.19E-01 


7.03E-03 


7.44E-06 


15 


4.00 


2400 


8.80 


524 


1.47E-05 


3.16E+01 


3.78E-01 


2.15E-03 


1.45E-05 


16 


4.00 


2400 


9.10 


524 


1.89E-05 


4.71E+01 


6.72E-01 


7.63E-03 


7.82E-06 


17 


3.85 


2600 


8.80 


393 


3.77E-06 


2.59E+01 


2.23E-01 


1.27E-03 


1.16E-05 


18 


3.85 


2600 


9.10 


393 


6.32E-06 


4.82E+01 


5.34E-01 


6.07E-03 


6.84E-06 


19 


4.00 


2600 


8.80 


524 


8.44E-06 


3.07E+01 


2.60E-01 


1.48E-03 


1.23E-05 


20 


4.00 


2600 


9.10 


524 


1.24E-05 


5.14E+01 


6.23E-01 


7.08E-03 


7.44E-06 


21 


4.00 


2800 


8.80 


524 


5.89E-06 


3.00E+01 


2.32E-01 


1.32E-03 


1.18E-05 


22 


4.00 


2800 


9.10 


524 


8.33E-06 


5.34E+01 


5.28E-01 


6.00E-03 


6.40E-06 


23 


4.00 


3000 


8.80 


524 


3.28E-06 


2.06E+01 


1.71E-01 


9.73E-04 


1.04E-05 


24 


4.00 


3000 


9.10 


524 


8.23E-06 


4.58E+01 


4.69E-01 


5.33E-03 


6.38E-06 



Table 5. Same as in Table[Tl but using the third {K-i) mean grain radius in the dust opacities. 



Mod. 


log(LJ 


Teff 


log(C-O) 


P 


<M> 


<"out> 


<fc> 


iPdlPg) 


<«gr> 




[Lo] 


[K] 




[days] 


[Moyr-'] 


[kms-'] 






[cm] 


1 


3.70 


2400 


8.20 


295 






4.73E-01 


6.43E-04 


2.85E-04 


2 


3.70 


2400 


8.50 


295 


2.54E-06 


1.57E+01 


8.18E-02 


2.22E-04 


1.03E-05 


3 


3.85 


2400 


8.20 


393 


- 


- 


5.15E-01 


7.00E-04 


2.96E-04 


4 


3.85 


2400 


8.50 


393 


5.18E-06 


1.73E+01 


7.37E-02 


2.00E-04 


1.18E-05 


5 


4.00 


2400 


8.20 


524 


- 


- 


6.27E-01 


8.52E-04 


3.73E-04 


6 


4.00 


2400 


8.50 


524 


9.03E-06 


2.19E+01 


8.25E-02 


2.24E-04 


7.65E-06 


7 


3.70 


2600 


8.50 


295 


5.75E-07 


1.43E+01 


7.14E-02 


1.94E-04 


1.78E-05 


8 


3.85 


2600 


8.20 


393 


- 


- 


4.99E-01 


6.78E-04 


2.49E-04 


9 


3.85 


2600 


8.50 


393 


2.53E-06 


1.74E+01 


5.74E-02 


1.56E-04 


6.73E-06 


10 


4.00 


2600 


8.20 


524 


1.91E-06 


3.40E+00 


1.86E-01 


2.53E-04 


2.01E-05 


11 


4.00 


2800 


8.20 


524 


- 


- 


- 


- 


- 


12 


3.70 


3000 


8.80 


295 


- 


- 


- 


- 


- 


13 


3.85 


2400 


8.80 


393 


5.84E-06 


2.90E+01 


2.48E-01 


1.41E-03 


1.22E-05 


14 


3.85 


2400 


9.10 


393 


1.09E-05 


4.30E+01 


5.74E-01 


6.52E-03 


7.20E-06 


15 


4.00 


2400 


8.80 


524 


1.44E-05 


2.94E+01 


2.93E-01 


1.67E-03 


1.24E-05 


16 


4.00 


2400 


9.10 


524 


2.01E-05 


4.71E+01 


6.60E-01 


7.50E-03 


7.93E-06 


17 


3.85 


2600 


8.80 


393 


3.61E-06 


2.90E+01 


2.71E-01 


1.54E-03 


1.12E-05 


18 


3.85 


2600 


9.10 


393 


6.57E-06 


4.76E+01 


5.25E-01 


5.96E-03 


6.79E-06 


19 


4.00 


2600 


8.80 


524 


8.68E-06 


3.13E+01 


2.40E-01 


1.37E-03 


1.18E-05 


20 


4.00 


2600 


9.10 


524 


1.33E-05 


5.06E+01 


6.14E-01 


6.97E-03 


7.45E-06 


21 


4.00 


2800 


8.80 


524 


5.31E-06 


3.01E+01 


1.85E-01 


1.05E-03 


1.04E-05 


22 


4.00 


2800 


9.10 


524 


8.10E-06 


5.54E+01 


5.23E-01 


5.94E-03 


6.25E-06 


23 


4.00 


3000 


8.80 


524 


3.58E-06 


1.98E+01 


1.44E-01 


8.20E-04 


9.75E-06 


24 


4.00 


3000 


9.10 


524 


8.12E-06 


4.91E+01 


4.68E-01 


5.32E-03 


6.32E-06 



